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Motivated by the classification ol non-equilibrium steady states suggested by R.K.P. Zia and 
B. Schmittmann in J. Stat. Mech. P07012 (2007), we propose to measure the violation ol the 
'detailed balance' criterion by the p-norm {\\K* 1 1 ) ol the matrix formed by the probability currents. 
Its asymptotic analysis, lor the totally asymmetric exclusion process, motivates the definition ol a 
'distance' from equilibrium K* obtained lor p — 1. In addition, we show that the latter quantity and 
the average activity {A") are both related to the probability distribution ol the entropy production. 
Finally, considering the open-ASEP and open-ZRP, we show that the current ol particles gives an 
exact measure ol the violation ol 'detailed balance'. 
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I. INTRODUCTION 



In nature, typically every system is governed by 
well-known 'deterministic' physical laws. However the 
microscopic details are usually unreachable and a full 
description of the system is impossible. The enormous 
number of degrees of freedom and the apparent 'chaotic' 
motion of the microscopic elements leads to rather 
complicated phenomenon. To face such situations, the 
physics of statistical mechanics has been built over 
the last two centuries. Our best approximation is to 
assume that the interactions between the microscopic 
elements occur according to some probabilistic rules. 
It follows that the natural reformulation of many-body 
problems takes the form of stochastic processes [1]. In 
addition, almost all systems are in interaction with the 
environment. It is only in particular limits (for defined 
time and spatial scales) that closed systems emerge. 
Generally, currents of particles, heat or magnetization 
are induced by the environment and testify to the 
non-equilibrium activity. If the presence of flux is a 
signature of a non-equilibrium steady state (NESS), the 
absence of macroscopic currents does not necessarily 
imply that the system is at equilibrium. By definition 
equilibrium is characterized, on a microscopic level, by 
the well-known 'detailed balance' criterion. The latter 
'balance' is breaking for any non-equilibrium state and 
imposes non-vanishing currents of probability flowing 
between connected states. 

If equilibrium systems are well described by the theory 
of ensembles, no global formalism exists for systems out 
of equilibrium. Recently, general relations for driven 
systems have attracted a lot of interest. One should 
mention the Kawasaki relation [3|, the Jarzynski and 
Crooks relations [3-6] and 'fluctuation relations' that 
include the fluctuation theorems of Evans- Sear les [3, Q 
and Gallavotti-Cohen d \l^. While analyzing NESS 
case-by-case, research works are traditionally focused on 
the probability distribution {P*} of the micro-states. 



This approach is usually relevant since almost all observ- 
ables can be extracted from there. However, the analysis 
of currents or average production of entropy requires the 
knowledge of the transition rates which complete the 
characterization of the steady state. Recently, Zia and 
Schmittmann [ll|, 113 suggested a general classification 
of NESS, where a complete description of the system 
is given by the distribution of the probabilities and 
probability currents {P*,K*}. This classification allows 
the identification of the transformations of the transition 
rates that leave the steady state invariant. Along these 
lines, a set of invariant quantities has been derived for a 
class of steady states driven by the boundaries dynamics 
|13| . From the probability currents a definition of the 
'Euclidean distance' from equilibrium has been proposed 
in [11| . In this paper we suggest to measure the violation 
of 'detailed balance' and define the 'distance' from equi- 
librium by the p-norm (||A'*||p) of the matrix formed by 
the probability currents. We give, for the periodic totally 
asymmetric simple exclusion process (TASEP), the exact 
expression of ||ir*||p. We show, in the thermodynamic 
limit, that the p-norm vanishes for any p y^ 1. This 
result motivates the definition of the 'distance' from 
equilibrium by ||iir*||j^ that we prove to be extensive for 
the periodic-TASEP, open-ASEP and open-ZRP. One 
should mention that the 1-norm was first defined in 
[lj| and used to measure the violation of the 'detailed 
balanced' criterion for different reaction-diffusion models. 

We start our paper with a presentation of the master 
equation which governs the time evolution of Markovian 
stochastic systems. Based on the analysis of the asymp- 
totic behavior of the p-norm ||if*||p, for the periodic- 
TASEP, we define the 'distance' from equilibrium. In 
section HI, we show that the average activity {A*) and 
the 'distance' from equilibrium K* are both related to the 
production of entropy. In section IV, for the open-ASEP 
and open-ZRP, we show that the current of particles gives 
an exact measure of the violation of 'detailed balance'. 
Finally our results are summarized in section V. 



II. GENERAL FRAMEWORK 



B. The p-norm as a 'distance' from equilibrium 



When describing Markovian stochastic many body sys- 
tems, the master equation is the most general equation 
that one has at his disposal. Its formalism is used to de- 
scribe in and out of equilibrium systems, such as chemical 
reactions, enzyme kinetics, biological populations etc... 
. In this section, we give a presentation of the master 
equation and introduce the matrix formed by the prob- 
ability currents. For the periodic- TASEP we discuss, in 
the thermodynamic limit, the asymptotic scaling of the 



p-norm 



\K* 



that motivates the definition of the 'dis- 



tance' from equilibrium by \\K* 



A. Master Equation 

Considering the continuous dynamics of a many-body 
system, we note X a particular state of the configuration 
space. The probability Px{t), to find the system in the 
state X at time i, is governed by the master equation 



dtPxit) 



Y=^X 



JW, 



with the net probability current 



K^{t)^w^xPY{t)-w^Px{t), 



(1) 



(2) 



where Wy is the rate of the transition from the configu- 
ration X to Y. We assume the dynamics to be ergodic 
which imposes that the system has a unique stationary 
state independent of the initial condition. The master 
equation ([T]) expresses the time derivative of the proba- 
bility as the balance of the currents flowing in and out 
the state X. Defining J+(X) and J-{X) as 

j+ix) = J2 <Py(.t)^ J-ix) = J2 ^yPxit), (3) 



Y^X 



Y^X 



the equation ([T} simply states the conservation of proba- 
bility through the expression dtPx{t) = J+{X) — J-{X). 
In the steady state, the current appears to be 'globally 
balanced' on each state so that J+{X) = J-{X). By def- 
inition, all equilibrium distributions verify the 'detailed 
balance' criterion 



i^*i=<P^-<Fj=0. 



(4) 



However, for NESS this criterion is locally violated and 
the probability current on each bond K*y is not zero. 
In the network, where nodes are states of the system and 
links are possible transitions between states, the violation 
of 'detailed balance' imposes the existence of loops of 
current. In such a situation the system presents generally 
(but not always) macroscopic currents. 



In the stationary state, a natural definition of the 'dis- 
tance' from equilibrium is given by the norm of the ma- 
trix defined by the elements K*^. For any systems evolv- 
ing in a finite configuration space, we define the p-norm 
by 



\K* 




p£ 



(5) 



where the factor 1/2 has been introduced for conve- 
nience. This quantity measures the violation of 'detailed 
balance' and vanishes at equilibrium only. In references 
[nl [l3 | the authors suggested to define an 'Euclidean 
distance' from equilibrium by ||ii'*[[2. In addition, they 
showed that each quadratic element |i4r*y p appears in 
the expression of the entropy production. 

Choosing the periodic- TASEP as an example, we ob- 
tain an exact expression of \\K* \\ as a function of the av- 
erage current of particles (j*)- Let us remind the reader 
that the latter model is defined on a one-dimensional pe- 
riodic lattice, where particles are strictly jumping for- 
ward to the next nearest neighboring site. The tran- 
sitions are only possible if the target site is empty and 
occur with a rate equal to one. Our work is facilitated by 
the fact that the stationary probabilities are equiproba- 
ble and that the dynamical process is totally asymmetric 
{w^ = if Wy = 1). Our calculations give for the p-norm 



\K* 






i/p 



(6) 



where fi is the dimension of the configuration space and 
L is the size of the system. As a consequence, while 
taking the thermodynamic limit and keeping the density 
of particles constant, we can show that ||-ft'*|L vanishes 
for all p ^ 1. Writing p the average density of particles, 
the asymptotic analysis [L » 1) gives us 



where 



|A'*|1„ oc Vii+I/P X e-'^^^i-i/p)^ 



-pln/3- (l-/3)ln(l ~ p). 



(7) 



(8) 



In this particular case, as in many others, the total cur- 
rent of particles {J*) appears as the simplest observable 
that characterized the non-equilibrium behavior of the 
system. This quantity increases with the system size 
while its density (j*) = {J*) / L remain constant. In 
our example, and with a natural understanding of the 
'distance' from equilibrium, the limit |lAr*|| — j- as 
L — >■ oo is clearly inadequate. The only 'satisfying' defi- 
nition consists of defining the 'distance' from equilibrium 
using the 1-norm. For simplicity, we will note hereafter 



K* = \\K* 



Exactly derived for the periodic- TASEP, 



open-ASEP and open-ZRP, K* appears to be extensive. 
Its relations to the activity and entropy production are 
presented in the following sections. 



III. ACTIVITY AND ENTROPY PRODUCTION 



If K*Y is well understood as the probability current, 
usually little is said about the role of the symmetric part 
A*Y = WyPx +w^Py and its properties. Introduced in 
[la lla], only recently has this quantity been analyzed, 
successively called trafhc fl7H19| , dynamical activity [20|- 
|23| or frenesy. Its central role in the non-equilibrium 
linear response theory and out-of-equilibrium dynamical 
fluctuation theory has been shown in [2J, [23 . Using the 
inequality | K*y \ < A*y , we trivially obtain the inequal- 
ity 



K* <{A*), 



(9) 



with (A*) = \ Ylx Y^x ^*Y which can be expressed as 
the average of the observable A{X) defined as 



A{X) = J2 



Wy- 



(10) 



Y^X 



The latter quantity corresponds to the escape rate and 
tells us how 'nervous' is the system in a given state. If 
one thinks about a system which presents, in its station- 
ary state, a macroscopic current (J*), the equation (|9|) 
should be understood as a generalization of the inequal- 
ity (J*) = ij;) + iJl) < \{Jl)\ + \{J1)\, where {j;) 
and {J'l ) are the averages of the currents flowing in op- 
posite directions. 

As already noticed in [ij, [I4I , the activity and the cur- 
rent of probability appear explicitly in the expression of 
total entropy production. Defining by {dS)r the average 
production of entropy measured over a time interval r 
(r ^ 1), the rate of creation a = {6S)r/T is given by 



E^- 



Wv\n 



XWy 



X,Y 



p 



X'^Y 

Py< 



(11) 



which can be explicitly written as 



1 . , / A*X 

i5:x*#in^^- 

x,Y 



K*^ 



A*§ 






(12) 



To pursue further we define the probability 'Pr(a, fc) to 
measure, in the stationary state, over a time interval r 
(t <IC 1), an activity a and a current of probability k. By 
definition, Vrici^k) is 



Vria, k)^22 PxWyT S {A*^ - a) (5 {K*^ - k) 

X,Y 



(13) 



Using the relations A*y — A*\ and K 
prove easily 



*x 

Y 



-K 



*Y 



{a + k)Vr{a,-k) = {a-k)Vr{a,k), 
which can be written as 
Vr{a,k) 



Vr{a,-k) 



s^^ with 5S = In (^) 
\a — kj 



(14) 



(15) 



The probability Vt{5S) of a creation of entropy 5S, over 
a time interval r (r <IC 1), is given by 



Vr{5S) 



+ OC 



da 



dk Vria, k) 6 



(55 -In 



a + k 
a — k 



(16) 



and such that equation p5)) leads to the well known 
ratio PriSS)/Pr{-6S) = e"^^. 

Particularly interesting results emerge from the analy- 
sis of the probabilities 



Pr{SS<0) = I d[5S)Vr{5S) 

J — CXD 

Vr{5S>Q) = / d{5S)Pr{5S), 



(17) 
(18) 



respectively, probability of a positive and negative cre- 
ation of entropy. From equation (J13l) one easily show 



PriSS^O) = ^ Pizw^r e (±ir" 



Y 1 



(19) 



X,Y 



where Q{x) is the Heaviside function defined by Q{x) — 
1 for X > and G(x) = otherwise. Expressed as a 
function of A*y and K*y , equation (|19p takes the form 

PriSS^O) = ^J2A*^[e{±K*^)+Q{TK*^)] 

X,Y 



X,Y 



(20) 



Restraining ourself to systems for which 'detailed bal- 
ance' is violated for every transitions such that K*y ^ 
yX ^ Y, one finally obtains 



PASS ^ 0) 



{A*)±K* 



(21) 



As a consequence, the probability Pr{SS = 0) is solely 
expressed as a function of the activity through the equal- 
ity 



Pr{SS = 0) = 1 - {A*)t. 



(22) 



In other words, the probability of a non-zero production 
of entropy is given by the activity produced over the time 
interval t. This is expressed under the equality 



(A*) 



PriSS^O) 



(23) 



Along the same line, the distance from equilibrium is 
measured by the difference 



K* = 



Vr{5S>Q)-Vr{5S <Q) 



(24) 



The latter equality can be explicitly expressed as the av- 
erage of the function sign 



K* = - I d{5S) Vr{5S)sign{5S), (25) 



where sign{x) is defined by 



-f a; < 

sign{x) = ^ a; = 

+ f a; > 0. 



(26) 



Since K* > we trivially obtain {sign{5S)) > 0, which 
simply states the inequality Vr[5S > 0) > VriSS < 0). 



IV. DRIVEN SYSTEMS 

If the calculation of the 'distance' K* can appear 
challenging, exact results are reachable for simple 
driven systems. In this section, we give the keys of the 
derivation which leads to the exact expression of K* 
for the open-ASEP and ZRP. These two systems have 
been the favorite toy models of the scientific community. 
They are some of the rare exactly solvable models, used 
as guides in the zoology of exotic behaviors found in the 
field of non-equilibrium many body systems. For these 
two cases, it is the particular structure of the stationary 
probabilities that allows an exact calculation of K*. 

As a preliminary result, from the inequality Q, one 
show that K* has to satisfy 



K* < \{j;)\ + \{ji)\, 

where (iZJ) are the average currents defined by 



(27) 



(28) 



X,Y 



where S ( X^Y 1 = 1 if the transition between the states 

X and Y is induced by the jump forward/backward of 
a particle. The total current of particles is given by the 
sum (J*) = (J!^) + (Jl). For such systems that do 
not present other reactions than the biased diffusion, the 
current of particles in the stationary state is directly re- 
sponsible of the production of entropy. Therefore a van- 
ishing current would imposes {5S) = so that we expect 
K* = K*{{J*)) with K*{{)) = 0. 



A. Open-ASEP 

The ASEP is the simplest one-dimensional driven 
model that includes the biased diffusion of hard-core par- 
ticles. Traditionally, the left and right diffusion rates are 
respectively denoted by q and p. Without loss of gener- 
ality, we are considering p > q which imposes (j*) > 0. 
The exclusion process imposes that the jump of a particle 
is only possible if the target site is empty. The system 
is driven at its boundaries by the interaction with two 
reservoirs of particles. On the left extremity, particles 
can be injected or removed with rates a, respectively 7. 
Identically, on the right boundary, the rates of injection 
and ejection are written S and /3. The open-TASEP is 
recovered for p — 1, q — and 7 = 5 = 0. It is well 
known that the exact expression of the stationary state 
can be expressed as the product of non-commuting matri- 
ces [2a,|23|. For a state X defined by the set of occupation 
numbers {771 (X), . . . ,'I]l{X)}, the stationary probability 
can be written as P 



X 



fL{X)/ZL with 



hiX) = (1^1 n ['7n(^)D + {1 - ^«(^)} E] \V), (29) 

ri=l 

and Zl = J2x fL{X). The matrices D and E satisfy 

pDE - gED = D + E = C, (30) 

and act on the vectors (W\ and \V) as 



(/3D-5E)|t/) = |F) 
(W^|(aE-7D) = {W\ 



(31) 
(32) 



It is important to note that the 'distance' from equilib- 
rium can be decomposed over the bulk {B*) and bound- 
aries contributions. We note K* = L* -\- B* -\- R* , where 
L* and R* are respectively the left and right boundary 
terms. Since for each transition that is induced by the 
jump forward of a particle [wy = p) the reversed tran- 
sition rate is w^ = g, the bulk contribution takes the 
form 



1 



^* = ^ E ip/(^) - '?/(^)i'^ (^^^ 



Zl 



(33) 



X,Y 



Considering a particle initially localized on the site n, the 
function f{X) and f{Y) can be written as 

/l(A) = (VK|$„(A)DE*„(A)|T/), (34) 

fL{Y) = (I^|$„(A)ED*„(A)|y), (35) 

where $„(X) = nr=iM'7»(^)D + {1 - ?7»(^)}E] and 
*n(X) = nt„+2 MX)Ti + {1 - 77,(A)} Ej. Using the 
commutation rule (j30p the contribution B* is given by 



B* = ^EE(^l*"(^)C*n(^)|V^)- (36) 



n X 



Since the bracket sandwich expresses a probability, each 
term of the sum is positive and the absohite value can be 
omitted. This last result leads to B* = {L- 1)Zl-i/Zl, 
where the ratio Z^^i/Zl corresponds to the average cur- 
rent of particles (7*). For details on the algebraic method 
see [26|,|27|. Along the same line, one has for each bound- 
ary, 

^* = ^El"/W-^/(^)l'^(^^^)' (37) 



with 



X,Y 



R* ^ ^Y.\Pf{X)-5f{Y)\5(xAY), (38) 



X,Y 



where 5 {X^Y\ , respectively 5 {X^Y\ , are equal to 1 

if the transition is induced by the injection of a particle 
on the left and the ejection of a particle on the right. 
Using the relations ((3T|) and ((32)) . one simply obtains 
L* = R* = (j*), such that 



K* = {L + l){f). 



(39) 



For p = q and for any set of parameters {a, /3, 6, 7} 
such that the stationary current vanishes, the system 
appears to be at equilibrium, characterized by the 'de- 
tailed balance' criterion. This situation is obtained when 
the density of the two reservoirs pa = a/{a + 7) and 
Pb = <5/((5 + /3) are equal {pa = pt = p)- In this situation, 
the steady state is described by the equilibrium (Gibbs) 
state of the lattice gaz at the density p [2g| . If the impli- 
cation K* = => (j*) =0 was obvious the latter relation 
states the equivalence K* = <^ (j*) = 0. 



B. Open-ZRP 

The zero-range process is a one-dimensional lattice 
model for which each site may be occupied by an ar- 
bitrary number m of particles. In the bulk, particles 
jump to the next nearest neighboring sites with rates 
9m function of the occupation number of the departure 
site. First introduced by Spitzer [20|, this model has re- 
gained interest since the observation of a condensation 
transition analogous to the Bose-Einstein condensation 
[131. See [ail-llll for a detailed study of the model. We 
are considering the open-ZRP where particles are added 
and removed through the boundaries. We note p9n and 
q6n the left and right transition rates associated to the 
move of a particle, with p > q such that (j*) > 0. On 
the first and last lattice sites, the injection rates of a 
particle from the reservoirs are respectively given by a 
and S. The ejection of a particle in the reservoirs, occur 
with rates q9n on the left boundary and p6n on the right 
boundary. The probability P^ of the state X, defined by 
{rii{X), . . . ,t]l{X)}, is factorized under the expression 



P*X^t[9niVn{X)), 



(40) 



lit, 

9nim) = ^Y[9^\ 



(41) 



fc=i 



where Z„ is the analogue of the grand-canonical partition 
function 



Zn^Y.-'n\[e-,'. 



(42) 



m k=l 



It was shown in ;34] that the fugacities z„ have to satisfy 
the equality 

pzn - qzn+i = a - 7Z1 = JSZL - (5 = {]*). (43) 

As done previously, we separate the bulk contribution 
[B*) to the boundary terms {L* and R*). Considering 
that the transition X ^^Y consists of the jump forward 
of a particle initially localized in n, we note m = 7y„(X) 
and m' = ?7„+i(X). The probabilities P'^ and Py take 
the form 

PI = 5„(m)5„+i(?7i') X W gkiVkiX)) (44) 

k^n,n-\-l 

Py = gn{m-l)gn+iim' + 1) X ]J gk{r]k{X)). 

k^n. 7i-\-l 

(45) 
In the expression of B* one has to evaluate the difference 



'.g„(m-l)5„+i(m' + l).(46) 



i.gn(™)g«+i(m') 



With the help of the relation ((43]) we derive the 
equality p9mgn{rn) = {{j*) + qzn+i)gn{m - 1) and 
q9m'+ign+i{m' + 1) = -{{j*) -pzn-)gn+i{m'). This fi- 
nally leads to B* — {L~ l)(j*) and, along the same lines, 
one show that the contribution of the two boundaries 
are both exactly equal to the current L* = R* ~ (j*)- 
As for the open-ASEP, the measure of the 'distance' 
from equilibrium appears to be given by expression p9p . 
One should mentioned that the periodic ZRP can be 
mapped to a periodic ID asymmetric exclusion process 
31, 32]. However, the new process under consideration 
is not as simple as the ASEP but defined by transition 
rates depending of the number of empty sites between 
particles and interpreted as long range interactions. For 
the periodic ZRP and ASEP one can prove K* = L{j*). 
Finally the factor L, for periodic systems, or L + 1, 
for open systems, is simply the number of possible 
transitions (between bonds and in/out of the system). 

Intuitively, for such one-dimensional systems, that do 
not present other reactions than the biased diffusion, we 
would conclude that a measure of the 'distance' from 
equilibrium and violation of 'detailed balance' is given by 
the calculation of the average macroscopic current. How- 
ever, a general derivation to any one-dimensional model 
remain a challenging problem since most of the time the 
stationary probabilities are unknown. 



V. OUTLOOK AND DISCUSSION 

Using the set {P*,K*} to classify the NESS, we pro- 
pose a definition of the violation of 'detailed balance' 
based on the p-norm of the matrix formed by the prob- 
ability currents. If, for finite system sizes, all norms are 
equivalent, we show by taking the thermodynamic limit 
of the periodic- TASEP that only ||ii'*||]^ appears to be 
a satisfying definition of the 'distance' from equilibrium 
{K*). Contrary to the entropy production, this quan- 
tity is invariant under the transformations proposed in 
[llL Il2l | that leave the NESS unchanged. In addition, 
we show that an upper bound of the 'distance' is given 
by the average of the observable A which measure how 
'nervous'/'active' the system is. Moreover, we show that 
both quantities are related to the probabilities of a pos- 
itive {Pr{5S > 0)) and negative {Vt{SS < 0)) creation 
of entropy SS measured over a time interval r << 1. 
This results lead explicitly to the expression of {A* ) and 
K* as a function of Vt{5S = 0) and {sign{5S)). For 
the open-ASEP and open-ZRP, the particular structure 
of the stationary state allows an analytical calculation 



which leads to the exact expression K* = [L + l){j*). 
Along the same line, more general results on the viola- 
tion of 'detailed balance' should be obtained in the near 
future, applying the matrix product technique to other 
systems that include charged particles or impurities. 
Even though the case-by-case analysis reveals K* as 
an extensive quantity, at this time, there is no general 
derivation of this result. We want to motivate the anal- 
ysis of the 'distance' K* , 'activity' {A*) together with 
the probability of entropy production Pt{SS) measured 
over an arbitrary time scale. We hope that this study 
will lead to a better understanding of the measure of the 
violation of 'detailed balance' and its direct relation with 
the entropy production. 
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